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a. 

C , Abstract 

C/3 , The purely classical counterpart of the Scattering Probability Matrix (SPM) j Sn,m f of the 

quantum scattering matrix S is defined for 2D quantum waveguides for an arbitrary number 
of propagating modes M. We compare the quantum and classical structures of | S„^rn P for 

C^ , a waveguide with generic Hamiltonian chaos. It is shown that even for a moderate number 

of channels, knowledge of the classical structure of the SPM allows us to predict the global 
structure of the quantum one and, hence, understand important quantum transport properties 
'■^ . of waveguides in terms of purely classical dynamics. It is also shown that the SPM, being an 

S ' intensity measure, can give additional dynamical information to that obtained by the Poincare 

O ' maps. 

Q ■ 



00 ■ 1 Introduction 

^^ . The S'-matrix is the most fundamental tool for analyzing quantum scattering phenomena in various 

f^ I fields of physics, for it provides us with the most complete scattering data M. Moreover, it is often 

^^ of interest to extend the analysis to the semiclassical regime. The first semiclassical formulations 

of the S'-matrix appeared in the early 70s in the works of Miller 0] and Marcus [pi in applications 

to atomic physics, and some years later extensions to their work were carried out by Heller B. 

C^ ' Their common approach uses the Feynman propagator in the WKB approximation, thus taking into 

account classical dynamics together with quantum mechanical interference, where the phases are 

given by the classical actions. In Ref. ||5[| some important further developments on the semiclassical 

treatment of scattering systems are listed. Quantum and semiclassical calculations of the 5-matrix 

Q ' have become essential for the understanding of transport phenomena in mesoscopic systems j^]. 

^ . In particular, in the ballistic regime, the conductance is well described by the Landauer-Buttiker 

formula G = ^J2nJ2m I ^n,™ P' where i„,m are the transmission elements of the S'-matrix. 

Semiclassical expressions for the transmission amplitudes for collinear leads were obtained by Jal- 

^ ' abert, Baranger, and Stone [||. See also Lin 
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One of the aims of studying ballistic motion in mesoscopic systems has been to relate the exper- 
imentally observed behavior of transport quantities, in the classical and quantum regimes, to their 
underlying classical dynamics ||l^ . This is particularly interesting when the associated classical dy- 
namics can be chaotic; then the purpose is to identify signatures of chaos in the transport pi[ . In a 
very recent example of this kind of work Ketzmerick |l2| showed that the fractal fluctuations of G as 
a function of magnetic field in a chaotic cavity are related to the Poincare-Birkhoff hierarchical struc- 
ture of the phase space of the corresponding classical motion. Previously, Jalabert, Baranger, and 
collaborators performed detailed quantum and semiclassical calculations of conductance in meso- 
scopic systems that display chaos in the classical regime pSJ. An important conclusion from their 
work is that the behavior of the average conductance can already discern whether the underlying 
classical dynamics is regular or chaotic. Specially relevant for our purposes here is their finding that 
the dominant contributions to the average quantum conductance are classical 113] . Thus it is natural 
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to expect that useful information may be obtained by analyzing purely classical quantities, disre- 
garding interference effects completely. Clearly, all the information contained in the phases of the 
quantum 5'-matrix, necessary to calculate e. g., the Wigner-Smith delay time, does not exist in the 
purely classical description. Nevertheless, as we shall show, important information can be extracted 
by studying the Scattering Probability Matrix (SPM); its elements are defined by the square modulo 
of the S'-matrix elements | Sn,m P, which give the transition probability for an incoming mode m to 
scatter into a mode n. The analysis of quantum and classical SP matrices is clearly relevant also for 
the study of the wave-ray correspondence of electromagnetic fields propagating in cavities |15|] since, 
under certain conditions, the wave equations are the same as for the quantum ballistic transport |16|| . 

In this paper we shall construct the purely classical counterpart of the quantum SPM valid for 
any 2D waveguide of arbitrary shape. Before doing so, in the next section we briefly review the 
definition of the S'-matrix in its application to cavities connected to leads. In Sect. Ill we construct 
the classical SPM and compare its quantum and classical structures for a model of a mesoscopic 
ballistic 2D wave-guide that displays generic chaos in the classical limit. We shall show that the good 
global correspondence between classical and quantum SPM enables us to understand the classical 
dynamical origin of features of the quantum SPM and to clearly identify the differences produced 
by the wave nature of the quantum state. In Sect. IV, we make some concluding remarks. 



2 ^-matrix for waveguides 

The S'-matrix relates incoming waves to outgoing waves, 

yout ^ S y™, (1) 

where F'" and 1/°"* stand for vectors specifying, respectively, waves coming into and going out of 
the interaction region. For a system composed of a 2D wave-guide of arbitrary shape connected to 
two leads, say left (L) and right (R) leads, the solutions in the leads are 

^^'^{x, y)^J2 K"" exp(zfc^'«x) + 6^-« exp(-zA:^-«a;)] cj^t^'iv), (2) 

m— 1 

where 



is the component of the wave function along the y-axis perpendicular to the direction of propagation 
(cc-axis) ; d^ stands for the constant width of the left lead which may be different from the width da 
of the right lead. For simplicity's sake we shall use dL = da for the rest of the paper. The sum is 
over all the propagating modes supported by the leads at a given Fermi Energy E. 

With this notation the S'-matrix and the incoming and outgoing waves can be written in the form 



'^ ' 1 X 7-jn I a \ -r TOUt ' 



The symbols t ,t', r, and r' in the S'-Matrix are M x M matrices, where M is the highest 
mode (given by the largest m beyond which the longitudinal wave vector kf^^ = ^/^ ^ ™j '^ be- 
comes complex). The symbols a^'^ and b^'^ stand for the vectors a^^ and b^^, m — I, 2...AI. The 
squared modulo element | tn,m P (| ^^,m P) gives the probability amplitude for a left (right)-incoming 
mode m to be transmitted to the right (left) lead into the mode n. Similarly, | r„^m P (| r'^ m P) i^ 



the probability for a left (riglit)-inconiing mode m to be reflected to the left (right) lead into mode n. 

The quantum SPM is simply defined as | Sn,m P; it gives the transition probability for the in- 
coming mode m to transmit or reflect into an outgoing mode n. 



3 Classical Scattering Probability Matrix (SPM) 

Since the energy of the system is given by its expression in the leads 

classically we can associate an angle Qm between the longitudinal component of the momentum hm 
and the total momentum ^2raeE/Ti. That is, 

rmrh 



sin ^ 



d^2mf.E 



(5) 



For a finite number M of modes there corresponds a range of angles A0m = 0.,„ — ^m-i for each 
mode m. The classical limit is M — oo. 

Consider a classical particle entering, say, from the left side and making an angle 9i, within a 
range corresponding to a given mode m. The particle (ray) will generally collide with the walls of 
the wave-guide a few times before exiting to the left or to the right, making a certain angle 9f^ to 
which we can associate a mode n ii 6f E Ad„- Initial conditions are specified not just by the angle 
but also by the initial position {x, y) along the left lead. In order to account for all possible types 
of trajectories, we take a large number of initial positions for each incoming angle Oi. By recording 
the number of particles scattered into the various ranges of 6 associated with different outgoing 
modes n, we obtain a distribution of outgoing modes for each incoming mode m. This distribution 
gives the classical counterpart of the matrix elements | i„.m P and | r„^m p of the quantum SPM. 
Similarly, to obtain the classical counterparts of | i^ „ p and | r^ „ p we repeat the above process 
but for particles entering from the right lead. This defines the procedure to construct the classical 
counterpart of the SPM. 

The wave-guide 

We now specify a wave-guide model on which to explore the quantum and classical structures of 
the SPM. We chose the geometry of the wave-guide to be that of a "rippled" billiard, shown in Fig. 
1, which is connected to two collinear leads of the same width. We have selected this particular shape 
for the wave-guide because it is known to display all features of an important class of dynamical 
systems, namely chaotic billiards that undergo the generic (Hamiltonian) transition to chaos |17|, [l^. 

Moreover, the finite version, depicted in Fig. 1, which serves as a model of a quantum or electro- 
magnetic waveguide, has been used to study certain transport manifestations of chaos in the classical 
pSJ as well as quantum |jl3| regimes. On the other hand, the infinitely long (i.e., periodic) version 
of the rippled billiard, introduced first in connection to beam acceleration problems pG], has been 
useful also for the understanding of typical features of crystals (e. g., energy band structure, LDOS, 
etc.) and their quantum-classical correspondence [Ell p2l p3| . 

Although, as a scattering system the finite version of the rippled channel is the relevant one, it is 
convenient first to review briefly the motion in the infinitely long rippled billiard, L — > cx). As usual, 
to get the dynamical panorama, we look at a Poincare map of the system. As the Poincare surface of 



section we choose, for reasons of symmetry, the bottom boundary y = 0; the Poincare map is given 
by the pair of Birkhoff variables {xj ^Oj), labehng the longitudinal components of the position and 
angle of the particle right after its j*'* collision with the bottom wall. Since the channel is periodic 
the Poincare map is on the cylinder {i.e., x is mod. 1). Depending on the geometrical parameters 
(average width 7 and amplitude u of the ripple) the dynamics is either regular, mixed, or fully 
chaotic. Figs. 2a and 2b show, respectively, typical Poincare sections for a wide (7 — 0.5, v — 0.12) 
and a narrow (7 = 0.25, v = 0.025) channel. In general, for small amplitudes of the ripple [v «1) 
wide channels (7 r^ 1/2) give rise to global chaos, whereas narrow channels yield mixed dynamics, as 
exemplified by Figs. 2a and 2b. For future reference we shall denote the system displaying globally 
chaotic dynamics (7 = 0.5, v — 0.12) as the G system and the mixed one (7 — 0.25, v — 0.025) as 
the M system. 

Since the Poincare plots of the periodic rippled billiard show topological chaos (i.e., a heteroclinic 
tangle) , it is not surprising that a finite rippled billiard connected to leads shows chaotic scattering, 
as evidenced by the fractal nature of its scattering functions, such as the dwelling time, reported in 
[ [l8[ . In fact, as is well known [|^, topological chaos is responsible for the fractality of the scattering 
functions. 



4 Results 

In the following we shall compare the quantum and classical SP matrices for both systems, G and 
M, and for various lengths of a rippled wave-guide. In all cases we will consider energies that allow 
for 33 propagating modes. 

The G-system 

Fig. 3a and Fig. 3b show, respectively, the quantum and classical SPM for the rippled wave-guide 
whose length equals one period of the ripple [L — 2n) . The resemblance between the classical and 
quantum SP matrices is remarkable. Let us consider first the reflection part of the SP matrices, 
say the left bottom block | r„_,„ p. Notice in the classical SPM the high intensity in the neigh- 
borhood of the {n,m) — (11,11) element and along a hyperbola like curve centered on it. The 
quantum SPM also shows the same pattern. The same is true for the cone-like shape starting at 
the {n,m) = (11,11) element. More impressive is the similitude of triangular shapes near the top 
right corner. This global correspondence enables us to predict, based solely on the classical pictures, 
important quantum transport features. For example the classical SPM predicts that there will be 
negligible reflection for modes m < 4. This is conflrmed in Fig. 4, which shows in detail how the 
incoming modes 2 and 3 do not reflect, whereas the incoming mode number 5 does reflect partially 
into the outgoing mode number 24, just as predicted classically [note the high intensity element 
(n, rn) = (24,5)]. These flgures also show that the modes 2, 3 and 5 transmit predominantly onto 
the same modes as the incoming ones; this would be just like the classical probabilities except that 
the quantum one shows, in addition, small transmission to some modes off the diagonal. Detailed 
analysis of the data shows that the classical SPM also gives transmission off the diagonal but it is 
not evident because their intensity is weak and almost uniform over all modes. This difference is 
due to quantum interference effects, which are responsible too for the larger width of the diagonal 
elements of the transmission parts. As another example, the classical SPM predicts that mode 11, 
incoming from the left, will reflect and transmit predominantly onto the same channel number, which 
is conflrmed by Fig. 4b. 

It is instructive to identify the type of trajectories that form the most salient features of the 
classical SPM since these are also evident in the quantum SPM. As an illustration, the triangular 
shape (see Fig. 5) that appears near the top right corner of the | r„_m p block results from incoming 
trajectories colliding only once with the rippled boundary in the neighborhood of a; = 1/2; the 



hyperbola-like curve and also the cone-like shape are formed by trajectories colliding twice with the 
rippled boundary. 

An important aspect of the quantum-classical correspondence which was not expected is the 
particle-like behavior that results from the interaction with the rippled wave-guide of certain plane 
waves. As an example, Fig. 6 illustrates this behavior for the incoming waves with mode numbers 
24 and 29. Note that a high intensity pattern is formed on the left side of the rippled wave-guide, 
resembling a ray trace. The angle this pattern makes with the horizontal, labeled a in the figure, 
corresponds precisely with the angle of reflection predicted by the classical SPM. In general we see 
that when there is a high intensity element in the classical SPM one can expect the wave function 
to form a ray pattern along the classical trajectory just outside the cavity. This may be regarded as 
a "short-lived scar" . 

In contrast, when the classical SPM shows homogeneous areas of low intensity probabilities, the 
quantum SPM is expected to show a mottled pattern of medium intensity probabilities. A homo- 
geneous area of low intensity classical probabilities results when incoming particles within a range 
A6m scatter uniformly throughout a much wider range of angles. This effect, the defocusing caused 
by the rippled boundary, is responsible for the strong sensitive dependence to initial conditions, the 
main ingredient of chaos. Clearly, the larger the number of periods forming the rippled wave-guide, 
the stronger this effect should be. Figs. 7a and 7b, showing the quantum and classical SP matrices 
for the same geometry as just above (7 — 0.5, v — 0.12) , but six times longer {L = 127r), confirms 
this expectation for the transmission parts. Comparison of the classical SP matrices (Figs. 3b and 
7b) shows that, with the exception of a high spot near the transmission element {n,m) = (28,28) 
and a short diagonal contribution (n, m) < (5, 5), all the distinguishing features of the transmission 
parts obtained for a one-period-long wave-guide (Fig. 3b) are washed out in the case of the six- 
period long wave-guide. The remaining high intensity diagonal elements for low modes are due to 
direct transmission, i.e., to trajectories that transmit without colliding with the upper wall. To get 
an estimate for the number of incoming modes that transmit predominantly onto the same mode, 
assume flat boundaries (since the amplitude of the ripple is small compared to the width 7) and 
consider a bundle of particles injected at {x,y) = (0.0,0.5). The particles can transmit directly (no 
collisions with upper or lower boundary) if their initial angle 9i is in the interval (—9c,0c), where 
9c = tan~^{2L/d), and L is the length of the channel. For the one-period-long wave-guide and six- 
period-long wave-guide, these angles are, respectively, 0.46 and 0.083 radians. Their ratio is 0.18, 
which agrees with the ratio between the lengths of the high intensity diagonals in the transmission 
of Figs 3b and 7b. 

We remark that the homogeneous spread of intensities in the transmission part of the classical 
SPM for (n, m) > (5, 5), is consistent with the "equal apriori distribution" of the S'-matrix required 
for the validity of the random S'-matrix theory approach to chaotic cavities pq | . But note that in 
the reflection parts of the classical SPM the inhomogeneity is especially strong. In fact, while the 
definite transmission structures of the one-period wave-guide have been somewhat washed out in 
the six-period wave-guide, classically, the reflection blocks remain practically the same as for the 
one-period-long wave-guide. This is because the reflection structures are mainly formed by particles 
reflecting within the first period of the ripple, see Fig. 5. On the other hand, backscattering after the 
first period of the rippled wave-guide shows up classically as an almost homogeneous spread of in- 
tensity throughout the reflection blocks (including the area below the hyperbola-like curve for which 
there was no reflection for the one-period wave-guide). Thus, while classically the definite pattern 
produced by the first period of the wave-guide persists, quantum interference of the backscattering 
from the whole wave-guide starts to destroy the pattern observed in the classical SPM. 

Another interesting feature comes from the analysis of the relatively bright spot observed in the 
transmission part of the classical SPM near the site (n, m) = (28, 28). A zoom of this spot is shown 
in Fig. 8a and a typical trajectory belonging to this pattern is shown in Fig. 8c. The distinctive 
feature of this type of trajectories is that they collide twice with the rippled boundary for each 



bounce with the flat boundary. These are periodic or quasiperiodic orbits advancing always to the 
right and form the stabihty island surrounding the stable period one fixed point shown in Fig. 8b. 
It is important to remark that this miniscule Poincare-Birkhoff structure (note the scale of the axis) 
is not visible in the whole Poincare map of Fig. 2a, even though its effect is clearly manifested in 
the classical SPM. Hence, we see that the SPM construction can give complementary information 
to that obtained by the Poincare maps because it is an intensity measure. On the other hand, such 
spot is not present in the quantum SPM because the size of the Poincare-Birkhoff structure is too 
small to be resolved quantum-mechanically (there are other spots visible but they do not correspond 
to the classical one; they originate from constructive interference). 

The M-system 

Now we examine briefly the classical and quantum SP matrices for the M system, Figs. 9a and 
9b show these matrices for a one-period-long wave-guide. Again, a quick comparison of these shows 
that the global features of the quantum SPM can be predicted by the classical SPM. We see that 
regions of high intensity areas in the quantum SPM correspond roughly to the areas of the classical 
SPM, albeit fluctuations within them. However, there are some important differences which we 
shall discuss now. Note that both classical and quantum SP matrices show that reflection occurs 
only for high modes but the classical reflection occurs only for modes higher or equal to to = 31 
while quantum reflection (although weak) exists even for modes as low as to, = 10. The mechanism 
responsible for the reflection of classical particles can be understood by examining the Poincare map 
of the infinitely long channel. Fig. 2b, which shows a large resonance island centered at a; = \. This 
resonance is produced by trajectories executing librational motion, bouncing between the two walls 
in the neighborhood of the widest part of the channel, x = j. It is clear then that particles entering 
the rippled wave-guide from the left at a; = can reflect (after one or several bounces) within the 
first period of the channel if their trajectories fall within the resonance island. Trajectories falling on 
the chaotic sea outside the resonance island (hence low transversal mode numbers) can also reflect 
via the chaotic separatrix but not within the flrst period of the ripple. The longitudinal momentum 
of these librational orbits is relatively small (see Fig. 2b), hence their transverse momentum is large. 
Detailed analysis using equation (5) and data from Fig. 2b shows that indeed the lowest mode that 
can reflect is to, = 31, in agreement with Fig. 9b. In contrast there is a strong quantum reflection 
for modes as low as to, = 25. Heisenberg's Uncertainty Principle is responsible for this difference, 
namely, the quantum state cannot resolve the flne classical boundaries defining the resonance island 
and consequently even smaller values of m can "tunnel" into the resonance island to cause partial 
reflection. 



5 Concluding Remarks 

We have studied quantum scattering properties of typical wave-guides with mixed and global chaos 
by examining the quantum Scattering Probability Matrix (SPM) and its classical counterpart. We 
emphasize that the definition of the classical SPM does not include any semiclassical aspects. We 
showed that the structure of the classical SPM allows us to predict the global structure of the quan- 
tum SPM. Since features of the classical SPM can be understood by analysis of the trajectories, 
it was possible to understand the classical dynamical origin of important features of the quantum 
SPM. Consequently, the analysis of the classical SPM of a given electron wave-guide system is useful 
for the understanding of its quantum transport properties, e. g., conductance. Plots of the classical 
SPM can be examined quickly to determine the influence of the cavity on the various modes. For a 
given energy, some modes may show ballistic behavior while others may display diffusive transmis- 
sion, as observed recently in Ref. |2(| . 

Our analysis of the quantum-classical correspondence of the SPM led us to discover the existence 
of "short-lived scars" . Specifically, we have seen that the wavefunction forms a ray pattern along the 



outgoing classical trajectory, for modes corresponding to high-intensity elements of the classical SPM. 
They are short-lived because after few bounces the ray pattern is destroyed by quantum interference. 

Finally, we wish to mention that certain small but relatively high intensity areas in the classical 
SPM lead us to discover the existence of extremely small Poincare-Birkhoff structures of the other- 
wise globally chaotic billiard. Hence, the SPM gives complementary information to that obtained 
solely by topological tools (e. g., Poincare maps) since it is an intensity measure. 
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FIGURE CAPTIONS 



Fig. 1 Geometry of the wave-guide. 

Fig. 2 Poincarc surface of section at 7 = for (a) {j,iy) ~ (0.5,0.12) and (b) (7,1^) = 
(0.25,0.025). 

Fig. 3 (a) Quantum and (b) Classical SPM, |5'„^mP, for the one-period wave-guide with 
(7,J^) = (0.5,0.12). 

Fig. 4 Reflection and transmission probabilities from the classical SPM for (a) the incoming 
modes 2, 3, 5, and (b) 11. The one-period wave-guide with (7, J^) = (0.5,0.12) is considered. 

Fig. 5 Typical types of trayectories that contribute to the zones marked in the reflection part 
k„,mP of the classical SPM of Fig. 3. 

Fig. 6 Wave function of the incoming mode (from the left) number (a) 24, and (b) 29. (a) 
a ^ 60°, and (b) a ^ 47° are the reflection angles predicted by the classical SPM of Fig. 3. 

Fig. 7 (a) Quantum and (b) Classical SPM, |S'„^mP, for the six-period wave-guide with (7, J^) — 
(0.5,0.12). 

Fig. 8 (a) Zoom on the transmission part of the classical SPM arround the site (n, m) = (28, 28). 
(b) Phase space generated by the trajectories that produce the structure of (a), (c) Typical orbit 
belonging to the pattern in (c). 

Fig. 9 (a) Quantum and (b) Classical SPM, |5'„^mP, for the one-period wave-guide with 
(7,1/) = (0.25,0.025). 
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